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We consider a complete list of simplified models in which Majorana dark matter particles anni¬ 
hilate at tree level to hh or hZ final states, and calculate the loop-induced elastic scattering cross 
section with nuclei in each case. Expressions for these annihilation and elastic scattering cross sec¬ 
tions are provided, and can be easily applied to a variety of UV complete models. We identify several 
phenomenologically viable scenarios, including dark matter that annihilates through the s-channel 
exchange of a spin-zero mediator or through the t-channel exchange of a fermion. Although the elas¬ 
tic scattering cross sections predicted in this class of models are generally quite small, XENONIT 
and LZ should be sensitive to significant regions of this parameter space. Models in which the dark 
matter annihilates to hh or hZ can also generate a gamma-ray signal that is compatible with the 
excess observed from the Galactic Center. 


I. INTRODUCTION 

Weakly interacting massive particles (WIMPs) repre¬ 
sent the most widely studied class of dark matter can¬ 
didates. WIMPs are motivated in large part by the fact 
that stable particles with weak-scale masses and inter¬ 
actions are generally predicted to freeze-out in the early 
universe with a density similar to the measured dark mat¬ 
ter abundance. Although neutralinos perhaps represent 
the most well-known example, viable WIMP candidates 
appear within a wide range of popular and well-motivated 
extensions of the Standard Model. 

In order for a stable particle species that was in thermal 
equilibrium in the early universe to avoid being produced 
in excess of the measured dark matter abundance, its 
number density must be efficiently depleted through self- 
annihilations (or through coannihilations). For particles 
in the GeV-TeV mass range, this typically requires an an¬ 
nihilation cross section of approximately uu ~ 2 x 10“^® 
cm^/s, as evaluated at the temperature of thermal freeze- 
out mm- In many cases, the Feynman diagrams respon¬ 
sible for dark matter annihilation can be related via a 
crossing symmetry to diagrams for dark matter elastic 
scattering with Standard Model particles. This has mo¬ 
tivated experiments designed to observe the elastic scat¬ 
tering of WIMPs with atomic nuclei. Over the past fif¬ 
teen years or so, the sensitivity of the direct detection 
experimental program has increased exponentially with 
time, strengthening constraints on average by a factor of 
two every 15 months. As a consequence of this progress, 
many otherwise viable dark matter candidates have been 
ruled out, as well as many classes of diagrams for dark 
matter annihilation. 

There are many scenarios, however, in which the dark 
matter can avoid being overproduced in the early uni¬ 
verse while yielding very low rates in direct detection 


experiments. Such possibilities include: 

• WIMPs that are depleted in the early universe 
through an efficient resonance mm- 

• WIMPs that are depleted primarily through coan¬ 
nihilations with another particle species [3 [S]. 

• Low-mass WIMPs (m^ ^5-10 GeV), which gener¬ 
ate signals below the detection thresholds of most 
direct detection experiments. 

• WIMPs with elastic scattering amplitudes that are 
kinematically suppressed by powers of velocity or 
momentum transfer. Dark matter mediated by a 
pseudoscalar is one well known example ISHS]. 

• Dark matter that annihilates to metastable parti¬ 
cles which subsequently decay to Standard Model 
states [MU]- In such scenarios, the dark matter’s 
elastic scattering cross section can be almost arbi¬ 
trarily suppressed. 

• WIMPs whose annihilations take place only 
through interactions with an odd number of Stan¬ 
dard Model states 

• Scenarios in which the abundance of dark matter is 
depleted as a result of a non-standard cosmological 
history, such as moduli domination and decay Eii- 
|28] or a period of late-time inflation [2M33 . 

• WIMPs that annihilate to Standard Model states 
which are not present in the nucleons that consti¬ 
tute the targets of direct detection experiments. 
In particular, the amplitudes for annihilations to 
leptons, heavy quarks, gauge bosons, and/or Higgs 
bosons lead to elastic scattering only at the loop- 
level. 
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It is the last of these possibilities that we consider 
in this paper, in particular focusing on Majorana dark 
matter that annihilates to hh or hZ final states. Al¬ 
though the elastic scattering of dark matter through 
loop-induced processes has been discussed extensively 
in the literature [9l [MH42] . and especially in the case 
of dark matter that annihilates to electroweak gauge 
bosons, such as wino-like neutralinos [551 - I57] . final states 
including Higgs bosons have been comparatively unex¬ 
plored. Such phenomenology can arise in a wide variety 
of theoretical frameworks, and dark matter candidates 
with these characteristics have become more attractive 
as direct detection constraints have become more strin¬ 
gent. Furthermore, such models can be motivated by the 
spectrum of the Galactic Center gamma-ray excess [la- 
I54| . which appears to be compatible with dark matter 
annihilating to hh or hZ final states [53] . 

In the following section, we describe the approach 
taken in this paper, which is based on a simplified model 
framework rather than any particular top-down theory. 
In Sec.|III[ we briefly discuss the Galactic Center gamma- 
ray excess and use this observation to motivate dark 
matter candidates that annihilates to hh or hZ final 
states. In Sec. 113 we describe the formalism for dark 
matter elastic scattering that we use throughout this pa¬ 
per. In Secs. [V] and VI we consider dark matter par¬ 
ticles that annihilate to hh (through the s-channel ex¬ 
change of a spin-zero mediator or the t-channel exchange 
of a fermion) or hZ (through the s-channel exchange of a 
spin-zero mediator, the s-channel exchange of a spin-one 
mediator, or the t-channel exchange of a fermion), respec¬ 
tively. We summarize our results and their implications 
in Sec. IVIII 


II. A SIMPLIFIED MODEL APPROACH 


In this paper, we do not adopt any particular ultravi¬ 
olet (UV) complete theory, but instead follow a bottom- 
up approach that makes use of simplified models. More 
specifically, we have assembled an exhaustive list of tree- 
level diagrams for the annihilation of Majorana fermion 
dark matter particles to hh or hZ final states, and cal¬ 
culated the annihilation cross section and the (loop- 
induced) elastic scattering cross section with nuclei in 
each case. This approach is similar to that taken in 
Ref. [7], in which we considered dark matter that an¬ 
nihilates to Standard Model fermions. In the case of 
annihilations to hh or hZ^ however, it is not as obvious 
that this could be the dominant annihilation channel in 
a UV complete theory, or that the same couplings re¬ 
sponsible for the annihilation cross section will provide 
the dominant contribution to the elastic scattering rate 
(see, for example, the discussion in Sec. 4.2 of Ref. [55]). 
In this section, we will discuss this issue and argue that 
such phenomenology could arise within the context of 
physically realistic models. 

We begin by considering a Majorana dark matter can¬ 


didate, X, which we take to be an electroweak singlet. 
The following effective field theory describes its interac¬ 
tions with the Higgs SU{2) doublet field, H: 

^ ^ *^p75)x 

+ -I- , (1) 


where A is a scale at which some electroweak-charged 
particle has been integrated out. After electroweak sym¬ 
metry breaking, the dark matter has the following inter¬ 
actions with the physical Higgs boson, h: 
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^A('^s+*Ap75)x/i-l-^y^x(A(,-l-iAp75)x/i^H- (2) 

If we were to assume that As ~ Ap ^ A), ~ A), ^ 0{1) 
and that < A^, the first term would yield large elas¬ 
tic scattering and annihilation cross sections through the 
CP-even and CP-odd vertices, respectively |7]. Such a 
large elastic scattering cross section would be in conflict 
with existing direct detection constraints. If, however, 
the couplings in Eq. [^were to respect the following hier¬ 
archy: 

As, Ap <C As, Ap X -y, (3) 


then the trilinear coupling between the dark matter and 
the Higgs will be subdominant, leading to rather different 
phenomenology. The leading interactions with the phys¬ 
ical Higgs particle will be a four-particle vertex in this 
case, leading to a low-velocity annihilation cross section 
(to hh) that scales as {av) ~ , and an elas¬ 

tic scattering cross section that is induced by a diagram 
involving a Higgs loop. 

In Eq.[2 one would also expect additional dimension-7 
operators involving the and Z, for example coupling 
the dark matter bilinear to and This 

would lead to terms in Eq. [^that couple xx fo ZZ, hZ 
and W~^W~. As these interactions are expected to be 
generated at the same scale. A, the effective coupling to 
hh could plausibly dominate without significant tunings 
or any large hierarchy between dimensionful scales. 

As an explicit example, we consider a simple super- 
symmetric scenario in which the dark matter is a ^100 
GeV bino-like neutralino, Xii and there exist two signif¬ 
icantly heavier, quasi-degenerate higgsino-like neutrali¬ 
nos, X 2 and Xa- Taking the large tan/3 limit, these col¬ 
lectively couple to the Higgs through the following in¬ 
teractions: £ D (v^gi/4)/iXiX2 + (v^5i/4)^Xi* 7®X3- 
As there is no Higgs-bino-bino coupling, tree-level elas¬ 
tic scattering is suppressed. The gauge eigenstates mix, 
however, leading to a non-zero Higgs—Xi — Xi coupling. 
After mass mixing, the relevant couplings are as follows: 


£ A 
+ 


{V2g,/i) /rx?X 2 + (^5i/4) 

gfvim^o - m^o ^ sin(2^)) 


2\/2m^o 


^XiXi- 


^2,3 


(4) 
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FIG. 1. The differential gamma-ray spectra for dark matter, x, annihilating to hh, Zh and bb with masses of = 130 GeV, 
110 GeV, and 65 GeV, respectively. It has been shown in the past that the Galactic Genter gamma-ray excess can be well-fit 
by dark matter annihilating to bb or hh for these values of the mass |53| . Given the similarity of the shapes (and not necessarily 
normalization) of these spectra, we expect that annihilations to Zh would provide a similarly good fit. 


For a modest hierarchy between the bino and higgsino 
masses, m^o ^ the hxiXi coupling can be effi¬ 

ciently suppressed (as can the ^XiXi coupling). In this 
case, annihilations will proceed dominantly to hh, and 
the loop-level diagram for elastic scattering will domi¬ 
nate over that from from the tree-level process (details 
of this calculation will be presented in Sec. VB): 


^loop 

<^tree 


1.3 X 10 


-4 


m^o 

A2 


0(3) 


m^o 

X2 

ITeV 
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In 


m^o 

X2 


100 GeV y 

J 


(5) 


Furthermore, it is well known that the supersymmetric 
parameter space contains significant regions, referred to 
as “blind spots”, in which the tree-level elastic scattering 
amplitude is suppressed by accidental cancellations [5B1 - 
[58] : loop-level processes could easily dominate in such 
regions. 

This example and others like it make it evident that 
UV complete and gauge invariant models can contain re¬ 
gions of parameter space in which dark matter annihila¬ 
tions proceed largely to pairs of Higgs bosons and elastic 
scattering is dominated by diagrams involving a Higgs 
loop. Furthermore, as direct detection experiments gain 
in sensitivity over the coming years, it is precisely these 
types of models (those with suppressed tree-level inter¬ 
actions) that will be the most difficult to test, and that 
will be the most interesting if no signal is detected. 

In light of these considerations, we proceed with an 
approach that is purely phenomenological in its philosos- 
phy. In the following sections, we explore simplified mod¬ 


els that introduce a minimal number of interactions, re¬ 
quiring only that they respect Lorentz invariance. Of 
course, UV complete models will generally introduce ad¬ 
ditional interactions, but since we parametrize our mod¬ 
els in a sufficiently general way, our results can be eas¬ 
ily mapped onto a given UV-complete scenario. In this 
sense, the cross sections calculated and presented here 
correspond to a phenomenological lower bound on the 
elastic scattering rate of WIMPs which annihilate di¬ 
rectly to Higgs bosons (in the absence of any strong can¬ 
cellations) . 


III. THE GALACTIC CENTER GAMMA-RAY 
EXCESS 


An excess of gamma rays from the region surround¬ 
ing the Galactic Center has been reported [13H53] . The 
morphology and spectrum of this excess are each in good 
agreement with that predicted from annihilating dark 
matter. More specifically, the best fits to the observed 
spectrum are found for dark matter that annihilates to 
quarks. For example. Ref. [53] find p-values of 0.35, 0.37 
and 0.22 for the fit of dark matter annihilating to bb, 
cc and qq, respectively. Dark matter annihilating to 
Higgs bosons also yields a fairly good fit to the data, 
with a p-values of 0.13 |53|. Annihilations to W'^W~ or 
ZZ provide a significantly poorer fit |53| (see, however. 
Ref. [55]). 

As far as we are aware, the Galactic Center gamma- 
ray excess has not been previously considered within the 
context of dark matter that annihilates to hZ. In Fig. [^ 
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we present the differential gamma-ray spectra from dark 
matter, y, annihilating to hh and Zh, with masses of 
= 130 GeV and 110 GeV, respectively. We note that 
for larger masses, the increased Lorentz boost between 
the dark matter and the hh ov Zh frames would result in 
significant broadening of the photon spectrum, worsening 
the fit to the Galactic Genter excess [IHHII]- In Fig. 
we present these spectra alongside that generated by di¬ 
rect annihilations to bottom quarks for = 65 GeV, 
which is known to give an adequate fit to the excess |53j . 
This figure illustrates that the spectral shape from the 
hh and hZ final states is very similar to that from direct 
annihilations to bb. Furthermore, the greatest departures 
are found at the highest energies, where the spectrum of 
the excess is least well measured. 


IV. DIRECT DETECTION FORMALISM 


For each simplified model considered in this paper, we 
will derive the effective Lagrangian describing the inter¬ 
actions of Majorana dark matter with quarks and gluons 
(following the formalism of, for example. Refs. pHl441) b 




u,d,s 

p(2) 




m: 




fll' 


X 




^(2) 


( 6 ) 


where is the scalar quark coupling, is the axial- 
vector quark coupling, is the scalar gluon coupling, 
and f!^'^ and are the quark and gluon couplings from 
twist-2 operators, respectively. The twist-2 operators are 
the symmetric traceless parts of the energy-momentum 
tensors and are given by: 


ilp^ q, (7) 

0% = 

In Eq. we have included all of the lowest-dimension 
terms that enable Majorana dark matter to scatter with 
nuclei without suppression by powers of momentum or 
velocity. Although the twist-2 gluon operator appears in 
this expression, its contribution to the effective nucleon 
coupling is subdominant compared to the scalar gluon 
coupling in all the cases that we consider throughout this 
study. We will neglect this contribution from here on. 

From the Wilson coefficients fq‘^\ and 

we can derive the spin-independent and spin-dependent 
scattering cross sections relevant for direct detection. 
The spin-independent cross section per nucleon is given 


by: 


0's/ 


/ m^rrin \ ^ 

lZfp + {A-Z)U 

\m^ +mnj 

A 


( 8 ) 


where ?7i„ is the nucleon mass and Z and A are the atomic 
number and atomic mass of the target nucleus, respec¬ 
tively. The effective couplings to protons and neutrons 
are given by: 


fp,Tl 

Tflpji 


E 


y(0) _|_ _Qip,n) j(2) 


q—u,d,s 



(9) 


where as is the strong coupling. For the mass fraction of 
light quarks and gluons, we adopt [60l [61] 0 

= 0.019, = 0.027, = 0.009, (10) 

= 0.013, f^'> = 0.04, 4”^ = 0.009. 

From these quantities, we can calculate the gluon con¬ 
tent, fxc'* = 1 ~ X) The quantities, 

q—u,d,s 

are the second moments of the proton and neutron parton 
distribution functions, which are scale-dependent. Since 
we will be performing our matching calculations at mh 
or mz, we will take the values as determined at these 
scales: 

mh : = 0.251, = 0.155, = 0.055, (11) 

mz : = 0.254, = 0.156, = 0.054. 

Values for the neutron are obtained by exchanging the up 
and down quarks in the proton case |64j . These quantities 
were derived using the parameterization and analysis of 
Ref. [Bij. 

The spin-dependent cross section of the dark matter 
with a proton or neutron is given by: 


12 / myTOn n 

crSD = — -—^— 

TT \m^ -I- rUp^n 


E 


q—u,d,s 


( 12 ) 


where the spin fractions of the proton and neutron are 
given by: = 0.77, = -0.49, A^f^ = -0.15, 

aL”’ = -0.49, = 0.77, and A^"^ = -0.15 [B5]. 

A complete treatment of these interactions would in¬ 
volve the following procedure. First, at the electroweak 
scale, one integrates out the Z, W^, t, and h by match¬ 
ing onto a complete basis of non-renormalizable operators 
(including ones that are kinematically suppressed and 


^ Other studies have favored significantly larger values for the nu¬ 
cleon’s strange quark content El E3]. The elastic scattering 
cross sections presented in this paper, however, are only impacted 
by this difference at the ~10% level. 
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FIG. 2. The diagrams for dark matter elastic scattering with nuclei corresponding to the case in which the dark matter 
annihilates to hh through the s-channel exchange of a spin-zero mediator, A (see Sec. V Al. 


that couple the dark matter to leptons). These opera¬ 
tors are then evolved down to the nuclear scale of roughly 
C>(GeV), including threshold corrections from integrating 
out the b and c quarks and the r lepton. Finally, once 
the operators are scaled down to ~ 1 GeV (at which 
the nucleon matrix elements are defined), one uses the 
nucleon matrix elements of the parton operators to gen¬ 
erate the effective nucleon couplings. For the simplified 
models considered in this study, however, the operators 
of Eq.j^are generated at leading order at the level of weak 
scale matching without strong cancellations, and thus the 
mixing of operators via the renormalization group flow is 
expected to only introduce a mild correction at the level 
of a few percent (see for example Ref. [35]). With this 
in mind, we derive the elastic scattering cross sections in 
this paper by performing a simple matching calculation 
at the weak scale (taken as mh or mz)- 


V. DARK MATTER ANNIHILATING TO hh 


of a spin-zero state. The Lagrangian for the simplified 
model describing this process is given as follows: 


£ D A X {>^8 + X + , (13) 


where A denotes the mediator and x is the (Majorana 
fermion) dark matter candidate. The annihilation cross 
section for this process is given by: 


crv = 


/4 

27r 



y(i-i:^)+A; 

(s - m\f + 


(14) 


where, rriA is the mass of the mediator and y/s is the 
center-of-mass energy. The width of the mediator, F^, 
potentially receives contributions from decays to hh and 
XX, given by: 


Majorana fermion dark matter can annihilate to a pair 
of Higgs bosons through two types of diagrams: the 
s-channel exchange of a spin-zero mediator, or the t- 
channel exchange of a fermionic mediator. In this sec¬ 
tion, we will calculate the cross sections for annihilation 
and elastic scattering in each of these scenarios. 


F(A^ M) = - 477^^)l/^ (15) 

(rrA — „ o ■> o t 

^ = 16nm\ ~ 


A. s-channel spin-zero mediator 

We begin with the case of dark matter that annihilates 
to a pair of Higgs bosons through the s-channel exchange 


Far from any resonances or particle production thresh¬ 
olds, the low-velocity annihilation cross section can be 
well-approximated by the first two terms of its Taylor 
series expansion, which in this case yields: 


av 



IGtt 



( 4 to ^ — m\Y 

Xp (l6m^ — 20m^m^ -|- m\mD — 2X1 

(4m2 — 


^a) {ml - ml) 


(16) 


Dark matter scattering with nuclei is dominated by the loop-diagrams shown in Fig.[^ These diagrams induce the 
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FIG. 3. For Majorana dark matter annihilating to hh through a spin-zero mediator, we plot the spin-independent elastic 
scattering cross section per nucleon as a function of the dark matter mass, for various values of the mediator mass, niA- In 
each case, the couplings, As,p, have been set to obtain the desired thermal relic abundance, = 0.12. In the left and right 

frames, we assume a purely scalar interaction (Ap = 0) and a mixed scalar-pseudoscalar interaction (As = Ap), respectively. We 
do not show results for As = 0, as the elastic scattering cross section is suppressed by powers of velocity in this case. The red 
regions in the upper portion of the frames are currently excluded by the LUX direct detection experiment , whereas in the 
blue regions we predict a cross section that is below the neutrino floor, making it difficult for dark matter to be detected by 
any planned direct detection experiment. Also shown are the regions within the projected reach of XENONIT (orange) and 
LZ (yellow) [67] . 


following effective couplings between the dark matter and Standard Model quarks and gluons (see Sec. IV): 


/r = 




IGTT^TO^m^ ’ 


(17) 


/r=i: 


asKvItJ-h 

32Tr^m'i 


dxidx2 ( “ l)a;ia:2 + ^^( 2^1 “ l)^xl{x2 - l)x2 




-m^(l - 2xi) {xi - l)xi{x 2 - 1) X 2 


where A = X 2 'm\ + {x 2 — 1)^^ and ^ = 3:2 + — 

1)(1 — X 2 )- Throughout this paper, j/g = —mqjv^ where 
V = 246 GeV is the vacuum expectation value of the 
Standard Model Higgs boson. 

In Fig. we plot the spin-independent elastic scatter¬ 
ing cross section per nucleon as a function of the dark 
matter mass, for three values of the mediator mass, niA- 
In each case, the couplings, Ag^p, have been set to ob¬ 
tain the desired thermal relic abundance, = 0.12. 

We compare these results to the current constraints from 
the LUX experiment |66| and the projected sensitivity of 
XENONIT (orange) and LUX-ZEPLIN (LZ, yellow) |S7]. 
The blue regions fall below the so-called “neutrino floor” 


(see, for example. Ref. (S^), where it unlikely that any 
direct detection experiment will be sensitive in the fore¬ 
seeable future. 

In Fig. 1^ we plot some of the phenomenological fea¬ 
tures of this model for the case of = 130 GeV and 
[ih = 300 GeV. For scenarios with a non-negligible value 
of Ap, the annihilation is dominantly s-wave (as seen in 
the first line of Eq. 16), which allows this model to gen¬ 
erate the observed normalization of the Galactic Genter 
gamma-ray excess (this criteria for this is taken to be 
that the low-velocity annihilation cross section is greater 
than 10“^® cm^/s |S3] and below the upper limit from 
Fermi’s observations of dwarf galaxies [68] . appropriately 
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rriA [GeV] tua [GeV] 



ruA [GeV] 


FIG. 4. For Majorana dark matter with a mass of 130 GeV, annihilating to hh through a spin-zero mediator, we plot the 
parameter space in which the predicted thermal relic abundance matches the observed dark matter density (solid black line) and 
the region in which the Galactic Center gamma-ray excess can be generated (green shaded band). The upper boundary of this 
region corresponds to the upper limit on the annihilation cross section from Fermi’s observations of dwarf galaxies [68] (green 
solid). The red regions in the upper-left portion of the first two frames are currently excluded by the LUX direct detection 
experiment [66], whereas in the blue regions we predict a cross section that is below the neutrino floor, making it difficult for 
dark matter to be detected by any planned direct detection experiment. Also shown are the regions within the projected reach 
of XENONIT (orange) and LZ (yellow) |67| . 

















FIG. 5. The diagrams for dark matter elastic scattering with nuclei corresponding to the case in which the dark matter 
annihilates to hh through the t-channel exchange of a fermonic mediator mediator, X 2 (see Sec. VBl. Crossed diagrams are 
not shown. 


rescaled for the case of annihilations to hh). In contrast, 
the annihilation cross section is velocity suppressed in 
the case of a purely scalar coupling (Ap = 0), leading 
to no appreciable indirect detection signals. In none of 
these cases are the current constraints from direct detec¬ 
tion experiments very restrictive, and unless tua is quite 
small, such scenarios predict elastic scattering cross sec¬ 
tions that are below the neutrino floor. Furthermore, if 
the coupling between the dark matter and the media¬ 
tor is purely pseudoscalar (As = 0), no detectable elastic 
scattering cross section is generated. 


B. f-channel fermionic mediator 

In this subsection, we consider Majorana dark matter, 
Xi, that annihilates to hh through the t-channel exchange 
of an additional fermion, X 2 ' 

C D h xi {K + X2- (18) 

The expanded low-velocity annihilation cross section in 
this case is given by: 


av ~ 


m2 \2\2 


2tt {-ml + + ml^) 


1 - 


mt 


m 


XI 


1 - 


1/2 


mt 


- 1/2 


967rm2^ (-m^-p m^^ + mlJ 
+ 6(Ap - A^) {ml - m^J + 3 - m^ J 

-8(A^-A^) m^i {ml-m\^f m\^ - {ml - m^^ 


{Xl-Xlfml^-{Xt-4XlXl + Xi)ml 


(Ap - 4A^A^ -h Xt) - 3ml^ {X^ - GA^A^ -h A^) 




3 (Ap - 4A^Ap -I- A^) ml + m^^ (“SA^ -I- ISA^Ap - SA^) 


+ 2(A^ - A^) m^i {ml - m^^f m^, + 3 {X^, + A^)^ m^^ (m^^ - ml) 


X2 


(19) 


Scattering with quarks and gluons occurs in this simplified model at leading order through the diagrams presented 
in Fig. Similar to Sec. |V A[ the full set of diagrams involves Higgs-vertex corrections, but now also includes the 
Higgs box diagram of Fig. [^b). Although this latter process is suppressed by the Yukawa couplings of the light 
valence quarks and is subdominant, for completeness we include it in our numerical analysis. The quark and gluon 
Wilson coefficients are given by: 


/r = 


-k 


yjrnq 

32'k'^v 

Vs I 


dx 3(a^-l)(A^ - A(+)m;^, (a: - 1)) 


IGtt^ 

2 


dX2 


Ai 

{1-X^-X2)il 

dxi - 


( 20 ) 


Ai 


2 1' 


A^~''^m^j (a;i(3a;2 — 4) — 2x2 -k 4) — 2A^ ''m^^{x2 — 2) 


-k 


m 


XI .^2 


A 2 


X 1 X 2 


X^ ^m^^ — X^~^'^m^.i^{xi — 1) 


f(2) = 

■><! 47^2 


rl—x2 


dX2 


dxi 


(1 - Xi - X 2 ) 


Ai 


A 


^^XIX2 


(A(- 


''X2 


- A(+)mx,(xi - 1)) 


-k -A^+^m;^iXiX2 
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4 °^ = E 


-a-sys 
3847r3z)mg Jq 


^ 3(a:-l)(A( - \^+'>m^^{x - 1)) 

Ai 


- l)a;i(a:2 - l)a;|(A(+)mxi + A( 

(^xl{x3 — 1) — xixz + a;i + X2X^ + \^~'’m^^{xi — l)a;i(a;3 — 1) 


+ ~ l)a^i(2;2 - l)a^3 

+ 2A2g4 ^^'^^"^xi(a^i - - l)x2{x^ - l)a;3 

+ 3A^3g2 - l)xi + 2)(x2 - 1)(X3 - 1)x3(A(+)toxi + 


+ 


^3S 

2 


3Aie3"“« 


ml{x2 - 1)(X3 - 1)X3 


A^+^TOj^i ^^3 [(xi - l)xi ((1 - 2xi)^ + 3 x 2 ) + 2 x 2 ] - (1 - 2xi)^(xi - l)xi^ 


+ A^ ^TO;^ 2 (xi - l)xi(l - 2 xi)^(x 3 - 1) 
1 


A3^4"‘Xi 


(xi - l)xi(X2 - 1 )X 2 (X 3 - l)a 


A*^+^m^^(xi - l)xiX2X3 + A*^ '’my^^m^^{xi - l)xiX2X^ 


- A(+^m^(l - 2 xi)^(x3 - 1) 


^A^+^TOy (xi - l)^Xi(X2 - 1)X2 (x 3 - l)x| 


Aie 


— 2xi)^(xi — l)xi(x2 — 1 )x2(x3 — lYx\{X ^'^"*+ A^ 






X^~''^m^img(l - 2 xiY(xi - l)xi(x2 - I)x2(x3 - 1)^X3^ . 


where we define the following: 


A(±) = A^ ± Xl, 

Ai = x(x - + : 


+ (1 - x)r 


A 2 = m^^(xi - l)xi +771^2^1 + ”2^(1 -xi - X 2 ), 


A 3 = 


_ 2:|x| 2 


^ + X 2 X 3 (m^j - + (x 2 - 1 )x 3 TO^ + (X 3 - l)mg, 


^ = xi(xi - 1)(1 - X 3 ) - X 3 . 


In the expression for above, the vertex diagram 
of Fig. I^a) contributes to the first line, whereas the box 
diagrams of Fig. [^b) contributes to the second and third 
lines. The first line of the expression for /g°^ corresponds 
to diagram [^c), while the rest of /g°^ results from dia¬ 
gram id). Although the box diagram in Fig. [^b) con¬ 
tributes negligibly to the overall scattering cross section, 
its 2-loop gluon equivalent, shown in Fig. [^d), has a 
significant effect due to the large top quark Yukawa cou¬ 
pling. 

In Fig. we plot the spin-independent elastic scatter¬ 
ing cross section per nucleon as a function of the dark 
matter mass, for two values of the mediator mass, rriy-^. 
In each case, the couplings, As,p, have been set to ob¬ 
tain the desired thermal relic abundance, = 0.12. 

Again, we compare these results to the current and pro¬ 
jected constraints from direct detection experiments. 


In Fig. we plot the annihilation and scattering be¬ 
havior for this model. As seen in Eq. |19[ s-wave anni¬ 
hilation can only occur if both the scalar (Ag) and pseu¬ 
doscalar (Ap) couplings are non-zero, and the annihila¬ 
tion rate is maximized in the case of Ag = Ap. In this 
case, as seen in the bottom-left and bottom-right frames 
of Fig. the required couplings are around 0(0.2), re¬ 
sulting in an elastic scattering cross section significantly 
below the neutrino floor. However, in the situation where 
either Ag or Ap vanishes, the annihilation amplitude is 
purely p-wave and necessitates 0(1) couplings in order 
to generate an acceptable thermal relic abundance. As a 
result, these larger couplings lead to scattering rates that 
might be within the reach of next generation experiments 
such as LZ. 

Among other features, the bottom-right frame of Fig. 
illustrates that special relations among the couplings Ag 
and Ap can lead to “blind-spots” in which there is an ac- 
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Xp=0 n^h2~0.12, As=0 




n^h^~0.12, Xs=Xp 



FIG. 6. For Majorana dark matter, xi interacting throngh a fermionic mediator, X 2 , we plot the spin-independent elastic 
scattering cross section per nucleon as a function of the dark matter mass, for two values of the mediator mass. In each case, 
the couplings, Xs,p, have been set to obtain the desired thermal relic abundance, = 0.12. In the upper left, upper right 

and lower frames, we assume a purely scalar interaction (Ap = 0), a purely pseudoscalar interaction (As = 0) and a mixed 
scalar-pseudoscalar interaction (As = Ap), respectively. The red regions in the upper portion of the frames are currently 
excluded by the LUX direct detection experiment [66], whereas in the blue regions we predict a cross section that is below the 
neutrino floor, making it difficult for dark matter to be detected by any planned direct detection experiment. Also shown are 
the regions within the projected reach of XENONIT (orange) and LZ (yellow) [67] . 
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m^j = 130 GeV, Ap=0 


m^j = 130 GeV, As=0 



Oh 





As 


FIG. 7. For Majorana dark matter with a mass of 130 GeV, annihilating to hh through a fermionic mediator, we plot the 
parameter space in which the predicted thermal relic abundance matches the observed dark matter density (solid black line) 
and the region in which the Galactic Genter gamma-ray excess can be generated (green shaded band). The upper boundary 
of this region corresponds to the upper limit on the annihilation cross section from Fermi’s observations of dwarf galaxies |68| 
(green solid). The red regions are currently excluded by the LUX direct detection experiment [66], whereas in the blue regions 
we predict a cross section that is below the neutrino floor, making it difficult for dark matter to be detected by any planned 
direct detection experiment. Also shown are the regions within the projected reach of XENONIT (orange) and LZ (yellow) | 67| . 


cidental cancelation between the contributions from the 
diagrams in Fig. Even for 0(1) values of the dark 


matter Yukawa couplings, this cancellation can persist if 
the combination of and Ap is tuned at the few percent 
level. 
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VI. DARK MATTER ANNIHILATING TO hZ 

Majorana fermion dark matter can annihilate to a Higgs boson and a Z boson through three types of tree-level 
diagrams: the s-channel exchange of a spin-zero mediator, the s-channel exchange of a spin-one mediator, or the 
t-channel exchange of a fermionic mediator. In this section, we will present the cross sections for annihilation and 
elastic scattering in each of these scenarios. 


A. s-channel spin-zero mediator 

Next, we turn our attention to Majorana fermion dark matter, y, that annihilates to hZ through the s-channel 
exchange of spin-zero mediator, A. The relevant interactions of this simplified model are described by the following 
Lagrangian: 


£ D H X (As -f x + QaZ^" {Ad^h - hd^,A) 

The annihilation cross section in this case is given by: 


( 21 ) 


9a 


IbTrm^ 


1 - 2 


ml+ml 




s — Idm^) + 5rn^ -|- s^ — lOsm 

( 22 ) 


where y/s is the center-of-mass energy. The width of the mediator, T^, potentially receives contributions from decays 
to hZ and xXj given by: 


T{A -A hZ) = 


r(A ^ xx) = 


9a 

167rm^ 


IbTrm^ 


m\ - 2m\ {ml + m|) -h {ml - to|) to| - 2 {m\ + ml) + 


1/2 


{m\-miy 


1/2, 


m\ - 4m^ J {m\ - 4m^) A^ -I- m^A^ 


(23) 


At low velocities and sufficiently far from resonance, this cross section can be expanded in powers of velocity: 
^ — 2(7m^-I-20m^)TO2-I-(3 to^-I- 4m^)^ 


(j^h - 2(to| + 4to2 )to^ -h (to| - 4m2 


(m\ — 4m^)2 




9a 


2567rm2TO^ (to^ — 4to^)^ 


2Apm^ 


8Apm^ {3ml - + 4m^) (ml - 2 (to| -I- 4m^) ml + (m| - 4m^)^ ^ 

— 2 {7ml + 20m^) m| -I- {3ml + 


212 
xJ 


{m\ — Ami) {m\ — 2 (m| -|- ) ml + {m\ — 4m^) 

X (m\{m\ — 2 (m| -I- 3m^) ml + m% + 8to^ — 6m|TO^) — 8m^ {ml — (2m| -I- 7 to^) ml + m% + 12m^ — 7m%m 

(ml - 2 (to| -h 4ml) + (w| “ 4^2)^ ^ (A^ -h A^) m^ (hm% “ 2 {7ml + 20 m 2 ) m| -F {3ml + 4ml)^^ |. 

In Fig. we show the dominant diagrams for elastic scattering in this case. This interaction is described by: 


-h 


£^ 


9Ays9q 

Jq 


dxi 


dx 2 ( 1 - 2 log ^ ) X (As -f Apiy®) x 1*7^9, 


(25) 


where A = ximl + X 2 m\ and = — g 2 T‘^/{2 cos Ow) is the axial piece of the Standard Model q-q-Z coupling. The 
7 ^ appearing in the quark part of this interaction guarantees that the resulting elastic scattering cross section will be 
both spin-dependent and suppressed by two or four powers of momentum for the As and Ap terms, respectively. As 
a consequence, the elastic scattering cross sections predicted in this model will be extremely small, and well below 
the neutrino floor. The reason that no qq term is present in Eq. 25 is that any interaction with two scalars and a 
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X 

q 



X 

q 


FIG. 8. The diagram for dark matter elastic scattering with nuclei corresponding to the case in which the dark matter 
annihilates to hZ through the s-channel exchange of a spin-zero mediator, A (see Sec. VIAI. 


Z is anti-hermitian, and thus cancels in the Lagrangian. As a result, only the scalar-pseudoscalar-Z vertex remains, 
providing a 75 in the effective quark bilinear. 

In Fig. 1^ we plot some of the phenomenological features of this model for the case of = 110 GeV and gA = 0.3. 
For scenarios with a non-negligible value of Ap, the annihilation is dominantly s-wave, allowing this model to generate 
the observed normalization of the Galactic Genter gamma-ray excess. In contrast, the annihilation cross section is 
velocity suppressed in the case of a purely scalar coupling (Ap = 0 ), leading to no appreciable indirect detection 
signals. 


B. s-channel spin-one mediator 


In this subsection, we consider Majorana dark matter, y, that annihilates to hZ through the s-channel exchange 
of a spin-one particle, Z': 




(26) 


The annihilation cross section in this case is given by: 


= -p-r-< (srriz' + 2rn^{m%i — Qsm\, + 

247rs3m|m2/— m^,)^-I-j i^ 

— 2 (s -I- to|) (^sm%, + 2771^(171%, — 6sm%, + 3s^)^m^ -I- s + I0sm% + s^) 

(m^/ — &sm\, + 3s^) m\ — 2s (idm^/ — 6 sm|/ + 3s^) m| -I- s^ {m%, — dsm^, + 3s^) 


(27) 


where y/s is the center-of-mass energy. The width of the mediator, F z ', potentially receives contributions from decays 
to hZ and XX'- 


T{Z' —>• hZ) — {^‘z’ + i'l^h ~ ■'^z) ~ 2m|, (m^ -|- m|) ^ 


1/2 


2_ (w|, -ml + to|)^ 




3/2 


247rTO2' 


Expanding the annihilation cross section in powers of velocity yields: 

,2N3/2 


av 




256Trm‘^m^m^ 


- (jnf^ — 2 (m| + 4m^) -I- (m| — 4 to^)" 


(28) 


(29) 


1/2 


9llJl{rni - 2 (to| + dm^) ml + (m| - dm^) 
30727rTO|TO^m2/ (w|/ — 4:mD^ 

— 7m^ -I- 20m^TO^ — 7m% + 32m^ -I- 2m| {7ml + 58m^) 


4 

mz, 







1 

+ 72 

mt-2 {ml + 2ml) 

2 1 4 A 2 2 

‘^h + '^z — 

ml,ml — 144 

ml — 2 {ml + 2ml) ml + m% — dm^m^ 

K 1 






J 
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rriA [GeV] 



tua [GeV] 



rriA [GeV] 


FIG. 9. For Majorana dark matter with a mass of 110 GeV, annihilating to hZ through a spin-zero mediator, we plot the 
parameter space in which the predicted thermal relic abundance matches the observed dark matter density (solid black line) 
and the region in which the Galactic Genter gamma-ray excess can be generated (green shaded band). The upper boundary 
of this region corresponds to the upper limit on the annihilation cross section from Fermi’s observations of dwarf galaxies |68| 
(green solid). In this case, we predict a cross section that is below the neutrino floor throughout the entire region shown. 


Scattering with nuclei occurs through the diagrams presented in Fig.[T0[ 
with these diagrams is given by: 


The dominant Wilson coefficient associated 


fW= ^dTlog I 

I67r2m|,n Jq x m\ -b (1 — a:) m\ Jq 



{xi -b X 2 + 1) 


(30) 


where A = xim\ -b X2m\ M 
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FIG. 10. The diagrams for dark matter elastic scattering with nuclei corresponding to th e case in which the dark matter 
annihilates to hZ through the s-channel exchange of a spin-one mediator, Z' (see Sec. VIBl. 


If we were to strictly adhere to our simplified model framework, we would use the Wilson coefficient given in Eq. 
to calculate the elastic scattering cross section in this model. In a realistic model with an interaction of the form 
however, mixing between the Z and Z' is an inevitable consequence of gauge invariance, leading to elastic 
scattering through tree-level Z exchange. For example, one might consider extending the Standard Model gauge group 
by a new broken 17(1), under which the Standard Model Higgs is charged. After electroweak symmetry breaking, the 
kinetic term for the Higgs leads to both an effective Z'-Z-h coupling as well as to mass mixing between the Z and 
Z'. As a result, these two effects cannot be decoupled, and a non-zero value for /i/j generically induces mass mixing, 
as well as an effective coupling between the dark matter and the Standard Model Z. This results in the following 
tree-level contribution to the axial-vector Wilson coefficient: 



where the mixing parameter, e, is given by: 


^■hV 

4(to|, - to|) ■ 


(31) 


(32) 


Mass mixing between the Z and Z' can also alter the value of the Z mass, which is strongly constrained by 
electroweak precision measurements. In particular, in the simple model described in the previous paragraph, the p 
parameter is shifted from its Standard Model value of unity as follows: 


p=l + e^ 



(33) 


Experimentally, p is bounded to be less than 1.0009 at the 95% conhdence level [TO], which translates into a fairly 
strong lower bound o n m^ /, even for modest values of ph- 

axial-vector quark couplings, fq^\ lead to spin-dependent scattering with nuclei. Although 


As discussed in Sec. 


IV 


spin-dependent scattering is comparatively weakly constrained, the cross sections predicted in this model are very 
large and well within the reach of current experiments. In the left frame of Fig. 11 we plot the spin-dependent elastic 
scattering cross section in this model, finding that it is strongly constrained by LUX |66] . In the right frame, we 
consider the case of = 110 GeV, and plot the parameter space in which the predicted thermal relic abundance 
matches the observed dark matter density and the region in which the Galactic Genter gamma-ray excess can be 
generated. The constraint from LUX excludes this simplified model, unless one considers dark matter that is very 
near resonance, 2m^ ~ mz'- We also include in the right frame the constraint from the measurement of the Z mass, 
as described in Eq. |33l 


C. t-channel fermionic mediator 

Lastly, we consider Majorana dark matter, yi, that annihilates to hZ (and necessarily to hh and ZZ) through the 
t-channel exchange of a fermion, X 2 '- 


C D {igl + 9x1^) X2 + hxi {Xs + Xpij^) X2- 


( 34 ) 
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ni^=110 GeV, //h=5 GeV 



FIG. 11. Phenomenology in the case of Majorana dark matter, y, annihilating to hZ through a spin-one mediator, Z'. In the 
left frame, we plot the dark matter’s spin-dependent elastic scattering cross section per nucleon as a function of the dark matter 
mass, for two values of the mediator mass. In each case, we have adopted a value of = 5 GeV and the coupling, g^, has been 
set to obtain the desired thermal relic abundance, = 0.12. In the right frame, we plot for the case of = 110 GeV the 

parameter space in which the predicted thermal relic abundance matches the observed dark matter density (solid black line) 
and the region in which the Galactic Center gamma-ray excess can be generated (green shaded band). The upper boundary 
of this region corresponds to the upper limit on the annihilation cross section from Fermi’s observations of dwarf galaxies |68| 
(green solid). The red region of this frame is currently excluded by the LUX direct detection experiment, when converted to 
apply to the case of spin-dependent scattering [66] . This constraint excludes the overwhelming majority of the parameter space 
for this simplified model. We also plot as a vertical dashed line the constraint from precision measurements of the Z mass, as 
described in the text. 


The expanded low-velocity annihilation cross section to hZ is given by: 

r 1 3/2 

m\ + (to| — — 2m^(4TO^^ -I- m|) 


CThZV 


1 

647rTO4^m| 


{ml - 2m2 - 2m2 -|- to|)2 


(m^i - m^^)Xsgl + {m^, + m^JXpgl 


1 2 


ml — 2 (m| -I- 4 to^J ml + (m| — 


(ml + m\ — — 2m2^) 


' (^hz{glD^ + Bhz{g^)^ + Chzg^g^ 


In the limit of m^^ 3> mh,mz- 


Ahz ^ TO® 


'X 2 '"-Xl 


I 4 4 


TD 4 4 

Bhz ^ 


+ 


Chz ^ 




1 — 4a; -b 3a;^ -I- 8a;® -b 3a;^ -b 4a;® -b a;® ) -b ^ml -b 


TO^^ ( 6 — 4a; -b lOa;^ -b 4a;^ ) -b 3 to^ -b 3 to| 




"XI 


"XI 


"XI 


6 -b 4x -b lOx^ -b 4x^ ) — 3ml — 3to| 

.4 - 5 , ^6 


1 -b 4x -b 3a;^ — 8x® -b 3a;^ — 4a;® -b x® ) -b ^ml + —m 


5 2 29 2 


X 2 


— 2 -b lOx^ — 2x‘^ — 6x® ) — ^ml — ^to| 


As Ap, 


( 35 ) 


where x = m^^ / m^^. 
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XI 


id) 



{h) 


XI 


FIG. 12. The diagrams for dark matter elastic scattering with nuclei corresponding to the case in which the dark matter 
annihilates to hZ through the t-channel exchange of a fermionic mediator, X 2 (see Sec. VICl. 


The velocity-expanded annihilation cross section to ZZ is given by: 


azzv 


47r 


2 \ 1/2 
z 


(’^X2 + ~ 


{«) 






247 rmT(-^7712(771^^ — i^x 2 ^xi ~ *^2) 

where in the limit of ^ mz- 


Snii m\ 8ml ml „ 

-- 

m% mi, 


X^zz + Bzz + Czz igx9x)^ ^ I 


(36) 


X2 ^Xi 
.6 ^6 


— ^Y2^yi 3 + 4x + + Sx^ + 7x'^ + 4x^ + 2x^ , 


Bzz ^ ^X 2 ^xi (^3 - 4x + 4x^ - + 7x^ - 4x^ + 2^*^ ), 

Czz ^ (j^xii^ ~ ~ ■ 


(37) 


Again, x = 771 ^^/ 777 ^ 3 . The annihilation cross section to hh in this model is the same as that given in Eq. 19 and we 
do not repeat this expression here. 

Turning now to elastic scattering, there are many diagrams to consider in this case, including the hh vertex correction 
(Fig. 12 7 , 12:), the hh box diagram |l^ , 12i), the ZZ vertex correction (Fig. 12:, 12 7 ), the ZZ box diagram 
[T^), and the hZ vertex correction (|l2^). The contribution from the hZ box diagram vanishes [10]. The resulting 
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Wilson coefficients are given by: 


f{ 0 ) ^ 1) - 2 tox 2 ) + -a;) - 2 to 


n ^ yj 

I Vs/mq 3 (x - l){X^-'>m^^ - X(+'>m^^{x - 1 )) 

+ 327r2r; J, Ai 


1 


L j•'^-X2 I 2 

dxo 


dx,{x,+x,-l)<^^ 


- 4 ^"1x2 (24 ^ + 9i'^Hx2 - 1)) 


+ {9^^\-^xxX2 + Xi + X 2 - 1) - 2g[ \xi - 


1 


^m\^m^^x\x2 + 2g[+'^m\^{xi - l)xlx2 


l>l-X2 




(1 - a:i - X 2 ) f 1 r 


Ai 


- (a^i( 3 a :2 — 4 ) — 2 x 2 + 4) — 2 A*^ ^'m^^{x 2 — 2 ) 


’^xi 2 


A*' — \^^'’m^^{xi — 1 ) 


( 38 ) 


5*'^^ 4 r^-x 2 f 3 1 r 

^ io ‘^2;i(xi +X2 - 1)< - ^4^^ + ^ \g^^'^m\^xi{xi + 1) - g^~'^m^^m^^xi 


Stt^v 


'“ r^-x 2 1 I r 1 r 

' dx2 I dxi —IXsg^ mx^{l-X2)+171x2 + Xpg^ rUxA^ - X 2 ) - ‘ 


'0 


'0 


^X 2 




2g^ '>171x1171x2x1x2 + 4g‘^'>mx^{xi - l)xlx2 
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— 1 W? n(+) ( _ 9 „(-) 


4°^ = X] Jo ^ - ^^(^1 - '^)xi{x2 - l)xl 4+^ {g 

^(xi - l)xi{x 2 - l)a ;3 2 g^^'>mx 2 Xi(^- xixo + a;i + 2:3 - 1 ^ 


X ^xi- 25x ’mx2 


.) 


+ (^a;i(a;3 - 1 ) - a;iX3 + Xi + ^ 2 X 3 ^ 

;(3(xi - l)xi + l)(x 2 - 1 )(X 3 - l)x^ (2g^g-'> 

xajCxi - 1 )xi(^4’^)(3(xi - l)xi + 1 ) - 2g^~\xi - l)xi) - 2g^~'^m^^ 


3A3^2-9\ 


2 




m (X2 - 1)(X3 - 1)X3 


( 3 ( 2^1 - 1)2^1 + 1) 2^2 ( 4 “^^ - 24 

(xi - 1 )xi(^4’^)(3(xi - l)xi + 1 ) - 2 g("){xi - l)xi) {g^'>m^^ - 2g^^m^^'j 




m^^{xi - l)xi(X 2 - 1)X2(X3 - l)x| 


4 '’g^x^'''m-lixi - l)xi(x3 - 1) +gq'^Jgx ^^Xl'^X 2 iXl - l)XiX2X3 


"^g^x ^rnl (3(xi - l)xi + (X 3 - 1 ) 


m^iTO^(xi - l)xi(X 2 - 1)X2(x3 - 1) X3 




gq ^( 2^1 - i)2;i - 24 ^ 1 x 1 ^ 2 ) 


(3(xi - l)xi + 1 ) - g[ 


5x^^W^iTO^(xi - l)xi(x2 - 1 )x2(x3 - 1 )^X 3 


Al^si'x "''xi'"'? 


gq \xi - l)xi - (^ 6 (xi - l)xi + 2 ^ ^ , 


where 


Ai = x(x — l)w^i + xm^^ + (1 - 
A 2 = xi(xi - + xito ^2 + 

A 3 = (Xi - l)xi + TO^^Xi + ; 

A 4 = + X 2 X 3 (m^^ - m 

As = X 2 (x 2 - + X2m‘^2 + 

^ = xi(xi - 1)(1 - X3) - X3 

4±) ^ (gir ± {gif 

5^^ - ^9qf ± {gff, 

± A^, 


and the vector and axial pieces of the Standard Model q- 
q-Z coupling are given by g^ = -g^ = g 2 Tjl/{2 cos 0w)■ 

In the expression for the scalar quark coupling, fjf\ 
the first and second lines correspond to ZZ and hh ver¬ 
tex corrections, respectively. Lines 3-5 and 6-7 of this 
expression result from the ZZ and hh box diagrams, re¬ 


x)to|, 

L - xi - X2)m|, 

-xi- X 2 ), 

J + (x 2 - l)x 3 m| -b (x 3 - l)m^, 
L — xi — X 2 )m'l + Xim|, 


spectively. The first and second lines of the expression for 
ff'^ arise from the ZZ box and hZ vertex correction, re¬ 
spectively. The quantities and ^ each result from 
the ZZ box diagram. 

In the upper frames of Fig. |13[ we plot the spin- 
independent elastic scattering cross section per nucleon 
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n^h2~0.12, Ap=5f"=0 


D.;^h^~0.12, As=Sf“=0 





[GeV] 


[GeV] 


FIG. 13. Phenomenology in the case of Majorana dark matter, Xi) annihilating to hZ throngh a fermionic mediator, X 2 - In the 
upper frames, we plot the dark matter’s spin-independent elastic scattering cross section per nucleon as a function of the dark 
matter mass, for two values of the mediator mass, and for either Xp = = 0 or Xs = g'^ = 0. In each case, we have adopted a 
value for the couplings to obtain the desired thermal relic abundance, = 0.12. In the lower frames, we plot for the case of 

nix = 110 GeV the parameter space in which the predicted thermal relic abundance matches the observed dark matter density 
(solid black line) and the region in which the Galactic Genter gamma-ray excess can be generated (green shaded band). The 
upper boundary of this region corresponds to the upper limit on the annihilation cross section from Fermi’s observations of 
dwarf galaxies |68) (green solid). The red regions are currently excluded by the LUX direct detection experiment |66j . while 
orange and yellow regions are within the projected reach of XENONIT and LZ, respectively |67| . Unless the dark matter is 
quite light, this model predicts cross sections that lie below the neutrino floor (blue regions), making it difficult for dark matter 
to be detected by any planned direct detection experiment. 
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as a function of the dark matter mass, for two values of 
the mediator mass, and for two selections of cou¬ 
plings: As = ^ 0 , Ap = g" = 0 and \ = g'^ ^ 0 , 

As = = 0 in the left and right frames, respectively. By 

setting these combinations of couplings to zero, we both 
minimize CP violation and approximately maximize the 
fraction of annihilations that proceed to hZ, as opposed 
to ZZ or hh final states. In particular, these combina¬ 
tions assure that the annihilation cross section to hh is 
p-wave suppressed, and that over much of the parameter 
space ahzvjozzv > ■ In each 

case, the non-zero couplings have been set to obtain the 
desired thermal relic abundance, ~ 0.12. While 

much of this parameter space predicts an elastic scatter¬ 
ing cross section that lies below the neutrino floor, this 
model could potentially lead to observable rates in future 
direct detection experiments if the dark matter is lighter 
than ~200 GeV (which is also required if the dark matter 
is to account for the Galactic Center gamma-ray excess). 

In the lower frames of Fig. we plot some of the 
phenomenological features of this model for the case 
of = 110 GeV (again for either Ap = = 0 or 

As = 5 “ = 0). For scenarios with a relatively light medi¬ 
ator, the annihilation cross section is dominantly s-wave, 
allowing this model to generate the observed normaliza¬ 
tion of the Galactic Center gamma-ray excess. Although 
current constraints do not yet restrict this model, future 
direct detection experiments are expected to be sensitive 
to this scenario. 


VII. DISCUSSION AND SUMMARY 

Direct detection experiments have become increasingly 
powerful in recent years. We are now able to test and 
exclude many models in which the dark matter couples 
at tree-level to Standard Model quarks. The null results 
of these experiments motivate us to consider dark matter 
candidates that instead couple primarily to gauge and/or 
Higgs bosons. Although such interactions can lead to ef¬ 
ficient tree-level annihilations in the early universe, they 
only generate a cross section for elastic scattering with 
nuclei at the loop-level. 

Past studies have considered the elastic scattering of 
dark matter through and Z loops in some detail. 
In this paper, we have instead focused on scenarios in 
which the dark matter is a Majorana fermion that cou¬ 
ples and annihilates to hh and/or hZ. In addition to 
generating small elastic scattering cross sections (con¬ 
sistent with existing direct detection constraints), dark 
matter that annihilates to hh or hZ final states can pro¬ 
vide an acceptable fit to the spectrum of the Galactic 
Center gamma-ray excess. Within a simplified model 
framework, we have assembled an exhaustive list of pos¬ 
sible descriptions for a Majorana dark matter particle 
that annihilates at tree level to hh or hZ. For each sim¬ 
plified model, we have calculated the annihilation cross 


section and the loop-induced elastic scattering cross sec¬ 
tion with nuclei. We have presented these cross sections 
in a format that can be straightforwardly applied to UV- 
complete models of interest. 

We have identified several particularly interesting sce¬ 
narios, including those in which the dark matter anni¬ 
hilates to hh or hZ through the s-channel exchange of 
a spin-zero mediator, or through the t-channel exchange 
of a fermion. Although we also considered spin-one me¬ 
diated annihilations to hZ, we found that such models 
are excluded by LUX unless the dark matter is very near 
resonance, m'^ ~ Future direct detection exper¬ 

iments such as XENONIT and LZ are expected to be 
sensitive to several of the cases considered here. In par¬ 
ticular, the fermion mediated models often yield cross 
sections that are within the anticipated reach of LZ. In 
contrast, models in which the dark matter annihilates to 
hh or hZ through a spin-zero mediator could plausibly 
exhibit elastic scattering cross sections that lie below the 
neutrino floor, making it possible that they will remain 
beyond the reach of direct detection for the foreseeable 
future. 

We note that several of the models considered here ex¬ 
hibit non-negligible low-velocity annihilation cross sec¬ 
tions, allowing for potentially observable indirect detec¬ 
tion signals (and for the ability to account for the Galac¬ 
tic Center gamma-ray excess). More specifically, each of 
the phenomenologically viable models considered in this 
paper yield significant low-velocity annihilation cross sec¬ 
tions, so long as the mediator’s coupling to the dark mat¬ 
ter has a non-negligible pseudoscalar component. Models 
with purely scalar couplings, however, can lead to p-wave 
annihilation amplitudes, and highly suppressed signals 
for indirect detection. 

In Table |Tj we summarize some of the most salient fea¬ 
tures of the models discussed in this paper. Given the 
breadth of the available parameter space, it is difficult to 
make simple statements that capture all of the possible 
phenomenology. Nonetheless, certain broad features of 
a given model can help us to project whether a future 
discovery with direct or indirect detection experiments is 
likely. 
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Final State 

Mediator 

Annihilation Cross Section 

Direct Detection Prospects 

hh (Sec. V) 

s-channel, spin-zero 

s-wave if As yf 0 

Difficult, but possible if As 7 ^ 0 
and rrimed < CJ(IOO) GeV 

t-channel, fermion 

s-wave if As 7 ^ 0 and Ap 7 ^ 0 

Possible if As Ap 

or As < Ap 

hZ (Sec. VI) 

s-channel, spin-zero 

s-wave if Ap 7 ^ 0 

Very Difficult 

s-channel, spin-one 

s-wave with comparable 
p-wave contribution 

Largely excluded by LUX in 
a UV-complete theory 

t-channel, fermion 

s-wave with comparable 
p-wave contribution 

Typically observable if 
moM ^ 200 GeV 


TABLE I. A brief summary of the prospects for indirect and direct detection in the case of Majorana dark matter that 
annihilates to hh or hZ. 
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